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GENERATION OF DIFFUSE NON-TBDEaiMAL ATMOSPHERIC PLASMAS 
Field of the Invendon 

This invention relates to the generation of gas plasmas, and more particularly to the 
15 electrically-eflacient production of diffuse non-theimal atmospheric gas plasmas. 

Review of the Art Known to the Applicant(s) 

Nonthermal gas discharges generated at atmospheric pressure find widespread use for 
ozone production (see B. Eliasson and U. Kogelschatz, 'TNfonequilibrium volume 

20 plasma chemical-processing", IEEE Trans. Plasma Science, vol. 19, pp. 1063 - 1077, 
1991), pollution control (see B. M. Penetxante, J. N. Barsley, and M. C. Hsaio, "Kinetic 
analysis of non-theimal plasmas used forpoUution control", Japan, J. Appl Phys. Vol. 
36, pp. 5007 - 5017, 1997.), and surfece modification of Polymer film g (see J. 
Friedrich, L. Wigant, W. Unger, A. Lippitz, and H. Wittrich, "Corona, spark and 

25 combined UV and ozone modification of polymer films WeBP23", Surface & Coatings 
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Technology, vol. 98, pp. 879 - 885. 1998.). Typically, such atmospheric plasmas 
consist of many nanosecond-scale microdischarges of streamer-like filaments (see A. 
C. Gentile and M. J. Kushner, "Reaction chemistry and optimization of plasma 
remediation of NxOy firom gas streams", J. Appl. Pkys. vol.78, pp. 2074 - 2085, 1995.) 
5 and much is known of their fundamental properties and flieir applications (see also M. 
Baeva, A. Dogan, J, Ehlbeck, A. Pott, and J. Uhlenbusch. "CARS diagnostic and 
modeling of a dielectiic barrier discharge". Plasma Chem. Plasma Processing, vol. 19, 
pp. 445 - 466, 1999 and R. Hackam and H. Akiyama. "Air pollution control by 
electrical discharges", IEEE Trans. Dielectrics and Electncal Instdation, vol.7, pp.654 
10 - 683, 2000.). More recently it has been suggested that nonfliermal atmospheric gas 
discharges can also be diffuse and luminous with the duration of their discharge current 
pulses in excess of sub-milUseconds (see E. E. Kunhardt, "Generation of large-volume, 
atmospheric-pressure, nonequihTirium plasmas", IEEE Trans. Plasma Sci., vol.28. 
pp.189 - 200. 2000; M. J. Shenton and G. C. Stevens. "Surface modification of 
15 polymer surfaces: atmospheric plasma versus vacuum plasma treatments", J. Physics. 
D: Appl. Phys., vol.34, pp. 2761 - 2768. 2001 and F. Massines. A. Rabehi. Ph 
Decomps, R. B. Gadri. P. Segur. and C. Mayoux, "Experimental and theoretical study 
of a glow discharge at atinospheric pressure controlled by dielectric barrier". Journal of 
Applied Physics, voL83, pp. 2950 - 2957. 1998.). Compared to their counterparts in the 
20 streamer-dommated mode, these diffuse nonthermal atmospheric plasmas have greater 
spatial uniformity, better temporal stability, and much lower gas temperature. typicaHy 
in the range 75»C - 150«C (see A. Schutze A, J. Y. Jeong, S. E. Babayan, J. Park, and 
G. S. Selwyn, "The atmospheric-pressure plasma jet A review and comparison to other 
plasma sources", IEEE Trans Plasma Sci, vol. 26, pp.1685 - 1694, 1998.). These 
25 properties make them particularly attractive for a number of key materials processing 



applications such as etching, deposition, and structural modification of polymeric 
surfeces. These diflEiise nonthermal atmospheric plasmas have been generated with 
dielecttically insulated electrodes at audio frequencies (see Massines et al supra), with 
un-insulated electrodes at radio frequencies (see Schutze et al, siq)ra\ and at 
microwave frequencies (see Shenton and Stevens, supra). Apart from their desirable 
spatial and temporal characteristics, diffuse nonthermal atmospheric plasmas have 
different electrical and chemical properties from that of their streamer-dominated 
counterparts. While many of then: ftmdamental properties remain to be fully 
understood, the focus of research has been largely experimental improvement of 
existing applications and empirical exploration of new applications. 

For controlled and optimised performance of material processing, it is always 
desirable to be able to adjust and tailor plasma properties. In the case of diffixse 
nonthermal atmospheric plasnaas, this has been pursued mainly through plasma rig 
deisigns and gas conqjosition, and to a lesser extent throu^ an appropriate choice of 
excitation frequency (See J. Park, 1 Henins, W. Hermann, and G. S. Sehvyn, "Gas 
breakdown in an atmospheric pressure radio-fiequency capacitive plasma source", J. 
Appl. Phys., vol 89 15 - 19, 2001). 

US Patent number US6228330 teaches the use of atmospheric-pressure plasmas for 
Ifae decontamination and sterilization of sensitive equipment and material, and provides 
evidence of their efficacy in inactivation of a number of bacterial spores. The patent 
describes a Systran for plasma recirculation, but discloses no shaping of the waveform 
of the applied voltage. 




The European patent application number EP 1 040 839 descibes a multi-stage 
process for sterilization using a gas plasma. The application discloses no shaping of the 
waveform of the applied voltage. 

US Patent number 6,118,218 describes an atmospheric pressure plasma tceater 
5 incorporating a porous metallic electrode. ITie application discloses no shaping of the 
waveform of the applied voltage. 

A number of features of non-thermal gas .plasmas are particularly important in 
regard to their applicability. First, the temperature of the gas plasma may be deleterious 
to objects with which the plasma comes into contact, for example during sterilization or 
10 surfece decontamination. Second, the power dissipation in the plasma has a direct 
bearing on the energy eflSciency of plasma generation, which is important for the 
overall economic efficiency of devices using the plasma, and is particularly important 
in applications remote from a readily-available supply of electrical power. 

Accordingly, it is an object of the present invention to provide a means to generate 
15 atmospheric-pressure plasmas at increased energy efficiency and at low temperature, 
whilst retaining the desirable feature of a high density of active species. 

Summary of the Invention 

In accordance with the inventive concept, these objects are addressed by the use of an 
20 applied voltage that exhibits a waveform (as defined herein) which is truncated (as 
herein defined) and/or which decays asymmetrically fix)m its peak value. 
Advantageously, the applied voltage, V, as a function of time, t, said time t being 
measured from any arbitrary instant, takes the form of a waveform, V(t), of cycle time 
T, wherein in at least one of the half cycles, i.e. . between (trn T) and (1?=i T + T/2) or 



between (t=i T + T/2) and (tF=(i+l)T), where i takes integer values, the wavefonn is 
characterised by the magnitude of the integral of the voltage with respect to time being 
greater in the first half of said half cycle than in the second half of said half cycle. 
Advantageously also, the applied voltage, V, as a function of time, t, said time t being 
measured from any arbitrary instant, takes the form of a waveform, V(t), of cycle time 
T, wherein in at least one of the half cycles, i.e. between (t=i T) and (t=n T + T/2) or 
between (iHi T .+ T/2) and (1r=(i+l)T), where i takes integer values, the waveform is 
characterised by a period of substantially constant voltage. 
Also advantageously, the applied voltage may be defined by equation El, below. 
Also advantageously, the applied voltage may be defined by equation E2, below. 
Also advantageously, the applied voltage may be defined by equation E3, below. 
Preferably, the applied voltage may be generated by the action of a control system, said 
control system using a measurement of the plasma discharge current as an ii^ut signal. 
Included within the scope of the invention is a method of generating a gas plasma 
substantially as described herein with reference to and as illustrated by the 
accompanying drawings. 

Brief description of the Drawings 

Figure 1 is a graph showing the reduced ionization coefficient of argon as computed 
with a Boltzmann solver (circles) and with the source term technique (+) together with 
reduced ionization coefficient of nitrogen using the source term technique (solid line). 
Figure 2 is a graph showing typical voltage-current characteristics of a helimn-nitrogen 
discharge under sinusoidal excitation with the gas voltage in solid curve, the discharge 
current in thick dashed curve, and the applied voltage in dot curve. 




Figure 3 is a graph showing a generalised peak4evelled wavefonru 
Figure 4 is a graph showing (a) A pulsed excitation voltage (solid curve) at lOkHz 
repetition firequency, made from a sinusoidal voltage (dashed curve) having its peaks 
levelled; and (b) time dependence of the peak-levelled excitation voltage (dot curve), 
5 the gas voltage (solid curve), and the discharge current (thick dashed curve). 

Figure 5 is graph showing normalized plasma power density (circles) and normalized 
electron density (diamonds) as a function of reduced magnitude of the applied voltage 
with 0.5% nitrogen. 

Figure 6 is a graph showing the normalized plasma power density (circles) and 
10 normalized electron density (diamonds) as a function of reduced magnitude of the 
appUed voltage with no nitrogen impurities. 

Figure 7 is a graph showing the normalized electron density (circles) and normalized 
metastable density (diamonds) as a function of reduced magnitude of the applied 
voltage with 0.5% nitrogen impurities. 
15 Figure 8 is a graph showing a generalised tail-trimmed waveform. 

Figure 9 is a graph showing (a) a peak-levelled excitation voltage at lOkHz repetition 
frequency (dashed curve) and with its pulse tail trimmed (solid curve); and (b) time 
d^endence of the peak-levelled and tail-trimmed excitation voltage (dot curve), the gas 
voltage (solid curve), and the discharge current (thick dashed cmrve). 

20 Figure 10 is a graph showing (a) a sinusoidal excitation voltage at lOkHz repetition 
frequency (dashed curve) and with its tail shaped with a Gaussian decay (solid curve); 
and (b) time dependence of the excitation voltage (dot curve), the gas voltage (solid 
curve), and the discharge current (thick dashed curve). 



The Inventive Concept 

The inventive concept considers pulsing plasma excitation voltage as a way to. 
control and improve properties of diffuse nonthermal atmospheric plasmas. For 
5 streamer-dominated atmospheric gas discharges, pulsed pi^ generation has been 
known to feciUtate better energy efficiency and greater control of the glow-to-arc 
transition althougji this has never been achieved for diffuse, non&ermal atmosphere 
plasmas. If similar improvements could be achieved for such plasmas, it would be 
usefbl for appUcations where plasma power consumption is an important issue, for 
10 example aircraft cloaking and industry-scale surfece treatment 

To illustrate the benefits of plasma pulsing in different plasma generation 
configurations, we will consider nonthemial almospheric gas discharges generated in a 
helium-nitrogen mixture and between two dielectrically insulated parallel-plate 
electrodes. Specifically the inventive concept attempts to answer (1) whether pulsed 
15 plasma generation can lead to saving in electrical power needed to sustain «ie generated 
plasma; and (2) if so how pulse shape and width may be adjusted to enhance energy 
efficiency. The approach uses on a one-dimensional computer code (See.X. T. Deng 
and M. G. Kong. "Parametric conditions for generation of stable atmospheric pressure 
nonlhennal plasmas", in the 2001 Annual R^ort, Conference on Electrical Imulation 
20 and Dielectric Phenomena 611 (ffiEE Catalogue Number: 01CH37225). pp. 677 - 680. 
2001) developed using a hydrodynamic model similar to those employed in most 
numerical studies for diffuse nonthermal atmospheric plasmas (See Massines et al, 
supra and F. Tochikubo. T Chiba, and T. Watanabe, "Structure of low-frequency 
heUum glow discharge at atmospheric pressure between paraflel plate dielectric 
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electrodes", Japanese Journal of Applied Physics, vol. 38, pp. 5244 - 5250, 1999). The 
key featuies of Ihe underlying model of difRise nonfliennal atmospheric plasmas as 
well as its numerical implementation will be described below, for clarity! Voltage- 
cuirent characteristics of nonthennal atmospheric plasma mider sinusoidal excitation 
5 will be used as an example to demonstrate Has utility of this approach. 

A physical model of nonOtemud atmospheric plasmas 

Difiuse nonthermal atmospheric plasmas may be induced and sustained between two 
parallel-plate electrodes, each one optionally coated with a dielectric layer and 
10 connected externally to a sinusoidal voltage source, typically at voltages in excess of 
IkV, and at audio frequencies. The. background gas is, for exanqjle, atmospheric 
heKum mixed with a smaU fraction of nitrogen (iq> to 1%) at room temperature of 
293K. and the dynamics of the generated nonthermal atmospheric plasma is described 
by the Boltzmann equation, well known in the art It has been established that the 
15 hydrodynamic assumptions can be applied to dif&ise nonlheimal atmospheric plasmas 
(See Massines et al and Tochikubo et al, supra). This reduces the Boltzmann equations 
to the continuity and momentum transfer equations for electrons and ions, thus 
faciKtating a macroscopic description of plasma dynamics without its microscopic 
details. On the other hand, the dynamics of electrons and ions are detemiined by the 
20 electric field in the space between flie two parallel electrodes, which has two 
components, one induced by the extemaUy appKed excitation voltage and the other by 
space charges. The electric field can be calculated by solving Poisson's equation. In the 
one-dimensional limit, fee Boltzmann equations and Poisson's equations are closely 
coiq>led together as follows 
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dt ^^^"^''^ Or Q^Tr (la) 

— ^=^K(^.0-«^(^,0] (Ic) 

where w+ and are the ion and electron densities respectively, and is the density of 
the neutral species considered. and iSi are the source terms for charged particles 
and neutral species respectively. jD and FT are difiRision coefficient and drift velocity 
with subscripts +, and i denoting respectively electrons, ions, and neutral species 
considered in the physical model. ^ is the electric field, \e\ is the charge of the electron, 
eo is the dielectric permittivity, and r is the variable representing the spatial position 
normal to the electrode. With the hydrodynamic assumptions, reaction rate coefficients, 
ionization coefficients, drift velocity, and diffusion coefficients can be approximated as 
a function of the electric field in the ionised gas between the two electrodes. Thus this 
physical model is similar to those used in Massines et al (stipra) and Tochikubo et al, 
(supra). The values of these coefficients are obtained firom relevant e>qperiments and 
several rate compilation studies which will be discussed in greater detail below. These 
allow eqs.(la) and (lb) to be solved to yield densities of electrons, ionic species, and 
metastables, which then allows Ihrou^ eq.(lc) calculation of electric field. Boundary 
conditions between the ionized gas and the dielectrically coated electrodes need to be 
carefully treated. To this end, a circuit equation is included to relate the electric field in 
the gas to the source voltage (the output voltage of the power supply) via a source 
resistor, i?s, and two serial capacitors each representing one dielectric coating layer. The 
numerical algorithm used to solve the above equations is essentially based on the 
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Patankar scheme, aod the discretization employs the upwind scheme (see S. Pataukar, 
"Numerical heat transfer and fluid flow ", Hemisphere publishing Co, 1980). 

The model* considers reactions involviag eight different species, namely (1) two 
ground-state neutral q)ecies, He(l^S) and N2(4^S); (2) two helium metastables, He(2^S) 
5 and He(2^S); (3) electrons; (4) ground-state atomic helium ions, He"*"; (5) ground-state 
molecular helium ions, 11^2 (2^SEu'^; and (6) ground state molecular nitrogen ions, 
Na'*'. In total, the model considers 17 reactions, as detailed in Table 1 below, together 
wift tiieir reaction mtes and reference sources. Products of these 17 reactions include 
two additional species, namely excited molecular heUum (He*) and atomic nitrogen 
10 (N)- Any subsequent reaction of these two additional species with the original eight 
species is not considered in the model. It is worth noting that the model is capable of 
evaluating production of molecular nitrogen ions that are not considered in reported 
numerical studies of difiRise nonthermal atmospheric discharges (e.g. Massines et al^ 
supra and Tochikubo et a/, supra). 

15 



Table 1: Reactions considered for a He-Nz discharge and their rate coefficients 



No 


Reaction 


Reaction rate 


Reference 


Direct ionization 


1 


He + e -> He+ 4- e + e 




1,2,3 . 


2 


N2 + e N2+ + e + e 




1,4 


Excitation 


3 


He + e -» He (23S) + e 




5 


4 


He + e He (2^S) + e 


A 


5 


5 


He (23$) + 2He (VS) ^ He2* + He (l^S) 


1.9Xl0^cmVs 


6 


De-exdtation 


6 


He (23S) + e He (l^S) + e 


4.2 XlO-^ cmVs 


7 


Penning ionisation 


7 1 


He (23S) +N2 He (l^S) +N2++ e 


8 X10-" cmVs 


5,8 


Stepwise ionisation 



) 
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No 


Reaction 


Reaction rate 


Reference 


8 


He (23S) + He (2^S) He (1 ^S) + He+ + e 


2.9 XlO^ ctnVs 


9 


Charge transfer 


9 


He+ + 2He (l^S) He2+ + He (l^S) 


6.3 X10^2 cmVs 


13,5 


Recombination 


10 


He+ + e He (VS) 


2 XI 0-^2 ^Vs 


10 


11 


Hea* + e Hea* 


5xlO-^*^cmVs 


11 


12 


He+ + e + He He (l^S) + He (2>S) 


1 X10*27 cmVs 


10 


13 


He2+ + 2e->He(l^S) + e 


7.1 Xl0-2° cmVs 


13.5 


14 


He2++ e -» 2He (23$) + He 


5 XI 0-^ cmVs 


6 


15 


Hea-^ + e +He He2' + He 


SxlO^^cmVs • 


11 


16 


N2+ + e N + N 


4.8 Xl0-« cmVs 


8 


17 


N2+ + 2e -» N2 + e 


1.4 XlO^^cmVs 


12 
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(3) J. Dutton, "A survey of electron swarm data", J. Phys Chem. Re£ Data, vdl.4, 
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helium afterglow at room temperature", Phys. Rev. A, vol. 13, pp. 1140 - 
1176, 1976. 

(8) A. Cenian. A. Chemukho, and V. Borodin, "Modeling of plasma-chemical 
5 reactions in gas mixture of C02 laser". Contribution to Plasma Physics, 

vol.35, pp. 273 - 296, 1995. 

(9) L.L.Alves, G.Gousset,andC.M,Ferreira,"Acomsional-radiativemodelfor 
microwave discharges in helium at low and intermediate pressures", J. Phys. 
D: Appl. Phys., vol. 25, pp. 1713 -1732, 1992 

10 (10) T. Quinteros, H. Gao, D. R. DeWitt, R Schuch, S. Pajek, S. Asp, and D. 

Belkic, "Recombination of D+ and He+ ions with low energy fiee electrons", 
Phys. Rev. A, vol.51, pp. 1340-1346, 1995. 
(1 1) P. C. Hill and P. R Hennan, "Reaction processes in a HeaYllu aaSg*) flash 
lamp", Phys. Rev. A., vol.47, pp. 4837 - 4844, 1993. 
15 (12) I. A. Kossyl, A. Yu Kostinsky, A. A. Matveyev, and V. P. Silakov, "Kinetic 
scheme of the non-equiUbrium discharge in nitrogen-oxygen mixtures". 
Plasma Source Science and Tech, vol.1, pp. 207 - 220, 1992. 
(13) Tochikubo etal, (supra). 

20 Ionization coefQcients in helium can be evaluated using Ward's fommla (A. L. Ward, 
x«^) or experimental data compiled by Dutton {supra). Both agree well with those used 
in Massines et al Qsupra), particularly at large reduced electric field (>50Td). Our 
model uses Ward's formula. On flie other hand, it is known from spectroscopic 
measurements ihat discharge dynamics in atmospheric heUum is very strongly affected 

25 by the presence of impuries, particularly nitrogen. For diffuse nonthermal atmospheric 
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helium dischaxges, Ifae presence of nitrogCTi impuries has been either ignored or 
implemented with its ionization coefficient approximated by that of argon. Since argon 
has a greater ionization coefficient that nitrogen, the approximation of using argon 
ionization coefficient overestimates- the ionization of nitrogen impurities. To this end, 
5 we note that it is possible to formulate analytically ionization source terms for 
discharges in argon, heUum, and nitrogen (See Peres, suprS^. Thus under the same 
conditions 



where a is Townsend's first ionization coefficient with its subscript A and B indicating 
10 respectively species A and B. If species A is argon and species B is helium, the above 
formula can be used to predict the ionization coefficient of argon from that of helium as 
calculated from Ward's formula. As shown in Figure 1, the argon ionization coefficient 
calculated with this technique, indicated by crosses 1 leads to an excellent agreement 
with that computed using a Boltzmann solver (see Massines et al, supra), indicated by 
15 circles 2. Therefore equation 2 offers a simple yet reliable way to estimate the 
ionization coefficient of one gas from the known ionization coefficient of a reference 
gas. Similarly if we assume species A is nitrogen and B is helium, the ionization 
coefficient of nitrogen is obtained and this is shown in Figure 1, indicated by the solid 
line 3, as a function of the reduced electric field. The ionization coefficient tiius 
20 calculated is used in our model for nitrogen (reaction #2 in Table 1). For excitation 
coefficients, we use that employed by Ben Gadri {supra). Similarly our model employs 
the same data used by Ben Gadri {supra) for diffiision coefficients and drift velocity. 
By way of example, and demonstration of the utility of the model presented above to 



_ ionization source term for species A 
ionization source term for species B 



(2) 
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accurately and correctly predict the evolution and properties of atmospheric diffuse 
nonthermal gas plasmas, a simulation is now presented, as follows. 

In this example, we consider a system with the two electrodes having a common radius 
5 of 2cm and separated with a gap fixed at O.Son. The total capacitance of the dielectric 
coatings on the two paiallel-plate electrodes is 70pF, and the source resistor in flie 
external circuit is 50a In this reference case, we consider a sinusoidal source voltage at 
lOkKfe with a peak voltage of 1.5kV. Before the power supply is switched onto flie 
plasma rig, the two electrodes are charged at 950V and so the initial gas voltage is 
10 950V. The background gas is 99.5% helium with 0.5%, nitrogen unless otherwise 
stated, and its breakdown voltage is assumed to be l.lkV, similar to the choice in 
Massines et al (supra). The secondary electron emission coefficient is dependent on the 
surfece condition of electrodes and as such it is not possible to choose a reUable 
coefficient that is appHcable to most cases. Numerically different secondary emission 
15 coefficients in the 0.01 - 0.2 range affect flie peak value of the discharge cunent For 
numerical studies considered here, we choose 0.2 for atomic helium ions, 0.1 for 
molecular helium ions, and 0.01 for nitrogen ions (see A. von Engel, "Ionized Gases", 
Chapter 3, reprinted by American Institute of Physics, 1994). Our numerical model 
assumes that diffuse nonthermal atmospheric plasma is established if its voltage and 
20 current are continuous and repetitive over at least 10 cycles of the applied voltage 
signal. Noting that discharge currents of both themial plasmas and streamer-dominated 
nonthermal atmospheric plasmas tend to exhibit microscopic current pulses in tens 
nanosecond or less, pulse duration of the discharge current is an additional indicator for 
diffuse atmospheric plasmas. Figure 2 shows a plot of the appUed voltage 4, the gas 
25 voltage 5, and the discharge current 6 as a function of time. Though not shown in 
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Figure 2, numerical resulte confirm that the discharge current 6 is in feet identical 
llirough many tens of cycles of the applied voltage 4 and so the discharge plasma is 
stable and repetitive tenqporally over a long period of time. The voltage-current 
characteristics are very similar to fliat obtained in a comparable experiment and its 
. 5 numerical simulation (Massines et al. supra), with a typical pattern of one discharge 
every half cycle of the appUed voltage. Both gas voltage 5 and discharge current 6 
exhibit very similar waveforms as those observed by Massines et al (supra) and 
Tochikubo et al {supra). The peak discharge current 6 in Figure 2 is between 52 - 58 
mA. lower than 90mA calculated by Massines etal.As shown in Figure 1 , this may be 
10 partly due to flie use of the larger ionization coefBcient of argon used by Massines et al 
to approximate the nitrogen ionization coefficient It is interesti^ to note that the peak 
current density of 7.2mA/cm^ (=90mA/«22) is the highest reported for difiuse 
atmospheric heHum discharges generated between dielectrically coated electrodes and 
at kilohertz firequencies. In a series of very similar e^qieriments performed at Tennessee 
15 University, the peak current density is found at most ImA/cm^ (J. R. Roth, "Industrial 
plasma engineering". Chapter 12, volA, loP publishing, Bristol, 1995.). A similar peak' 
current density of about 0.4mA/cm* (=11.25mA/3r32) was measured with 3kHz plasma 
excitation by the OkazaM group at Sophia University in Japan (T. Yokoyama, M. 
Kogoma, T Moriwaki, and S. Okazaki, 'The mechanism of the sterilisation of glow 
20 plasmas at atmospheric pressure", J. Phys. D: Appl. Phys., vol. 23, pp. 1125 - 1128, 
1990), whereas another Japanese group measured l.SmA/cm^ with a lOOkHz power 
supply. More recently a heUum discharge experiment performed at Minnesota 
University recorded a peak current density measured at O.SmA/cm^ with 15kHz driving 
voltage (L. Mangolini, K. Orlov, U. Kortshagen. H. Heberlein, and U. Kogelschatz, 
25 "Radial structure of low-frequency atmospheric-pressure glow discharge in heHum", 
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Appl Phys. Lett., vol. 80, pp. 1722 - 1724. 2002). Our computed cuirent density of 
4.4mA/cm2 (=55mA/;r2^) falls between Ifaese e:q)erimental measurcmenis, indicating a 
good cunent prediction capability of our model. 

Electric power consumption in the plasma is found to be 298mWW from our , 
5 model. TTiis is almost identical to SOOmW/cm^ measured by Massines et al ^supra) and 
similar to nimSNIct^ measured by Chen a/ (Z. Chen, J. E. Morrison. R. Ben Gadri. 
and J. R. Roth, "A low-fiequency impedance matching circuit for a one atmospheric 
uniform glow discharge plasma reactor", paper 6P63. presented at the 25* IEEE 
Jntemational Conference on Plasma Science, Raleigh, USA. June 1998). We have also 
10 calculated electron density and ion density across the space between the two electrodes. 
The time and spatial variation of electron density is similar to fliat computed in 
Massines et al (supra), A rapid rise of electron density occurs near 250;is and is a result 
of fee gas ionisation (or discharge) in that half cycle (250 - 300//s). and the relatively 
• uncharged electron density peak during the remaining period of the half cycle suggests 
15 a phase of inactive electron production. The peak electron density predicted with our 
model is 0.95x10'° cm\ not very dissimilar to 2xl0'° cm"^ by Tochikubo et al (supra) 
but much lower than-SOxlO'" cm-^calculated by Massines et al (supra). Again the large 
electron density reported by Massmes et al (supra) may be a result of their 
overestimated ionization coefBcient for nitrogen (see Figure 1). For difiiise atmospheric 
20 helium plasmas, there is very Uttle data of direct measurement of electron density for 
comparison with our calculated electron density. However indirectly the ion-trapping 
mechanism proposed by Roth (supra) may be used to estimate the average electron 
density and this can be captured by the following equation: 



n. -■ 



c* av 



2»2 (3) 



") 
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where n« is the average electron density* m electron mass, e electron charge, Vc electron 
collision frequency in helium^ Pov the average electric power density consumed in the 
plasma, and £b the peak electric j&eld in the plasma. Assuming Pav =300mW/cm^, 
jEb=3kV/cm, and Vc =1.8x10^^ Hz (see Roth, supra), eq.(3) yields 4.3xl0*cm"^. Given 
5 that the peak electron density is approximately 10 — 50 times greater than the calculated 
average electron density, the above calculation suggests that the peak electron density is 
between 4xl0^cm"^ and 2xl0^°cm"^. Thus our calculated peak electron density appears 
reasonable. This is also confirmed by electron density estimate using the current density 
formula of 7 - enefiE (where /i is the electron mobility) as employed for measuring 

10 electron density in micro-hollow cathode discharges (see R. H. Stark and K. H. 
Schoenbach, *T)irect current high-pressure glow discharges". Journal of Applied 
Physics, vol. 85, pp. 2075 — 2080, 1999). In graeral our simulated voltage-current 
characteristics, the peak dischiarge current density, the plasma power density, and the 
electron density yield a favourable comparison with available experimental and 

15 numerical data of comparable diffuse atmospheric heliimi discharges. While there is 
much scope to include other features of dififixse atmospheric plasmas, for example 
multi-dimensional effects, our plasma model described above is capable of capturing 
the main plasma features accurately and as such it will be used to explore tiie benefits 
of plasma pulsing. 

20 

The underl3dng inventive concept, realising a benefit of energy saving is best 
understood firom the voltage-current characteristic m Figure 2 (Massines et al, supra). It 
is seen that the discharge current 6 is of very low magnitude aft^ its peak between 
207//S - 245;zs in the half-cycle from 200;/s to 250//S, when the applied voltage 4 is 
25 very la/ge. This suggests that between the instant of this discharge cmrent pulse and 
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the polarity reversal point (that occurs at around 250jus) of the appKed voltage 4 in the 
same half cycle the applied voltage 4 does not necessarily contribute to plasma 
generation but may adversely heat up electrons. Through energy transfer ficom electrons 
to heavy particles, this undesirable electron heating would increase gas temperature. 
Therefore by reducing the applied voltage 4 over flris "inactive" period of electron 
production, flie gas temperature may be lowered and equally importantly the input 
power can be reduced. 

Description of the Preferred Ernhndimentfi 

We will now describe preferred embodiments of the invention by means of four 
examples. The first uses a waveform essentially comprising a peak-levelled sinusoid. 
The second employs a 'tail-trimmed' sinusoid, and the third employs a *tail-shaped' 
sinusoid. The fourth example describes how a control system can be used to shape the 
waveform of the applied voltage in response to the resultant discharge current 

Example 1-PeakrLevelled Siniisoid 

We first consider a simple pulsed excitation voltage waveform, shown in Figure 
4a, which is essentially derived from a sinusoidal wave with its peak (positive and 
negative) levelled to a flat top. Voltage waveform of Figure 4a can be easily obtained 
electronically, and so it represents a realistic option. By reference to Figure 3, the 
general equation describing the peak-levelled pulse waveform, Vft), is as follows: 
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F(/)= 



where : 



sm(2;r / 1) elsewhere 



2 



and: 

rem(2 7C f t, 2n) gives the isemainder after the division (2 ^ / {)l2a . 

5 

With reference to Figure 4, the magnitude of Ae original sinusoidal voltage, V„ is fixed 
at 1.5kV, whereas that of the peak-levelled voltage 7, varies from 0.4kV to l.SkV. 
The repetition frequency of the excitation voltage remains at lOkHQz. 

Figure 4b shows the voltage and cunent characteristics of a diffuse nonthermal 
10 atmospheric gas discharge generated with = 1 .2kV. The discharge current 8 remains 
* repetitive and stable through many tens of cycles of the applied voltage 7. Therefore 
under the conditions considered here, pulsing the plasma-generating voltage 7 does not 
appear to significantly affect the estabhshment of diffuse nonthermal atmospheric 
plasma. More specificaUy the waveform of the discharge current 8 remains relatively 
15 unchanged from that under the sinusoidal excitation, and the peak current is again 
around 55 mA very similar to the sinusoidal case of Figure 2. 

The basic volfage-current characteristics remain relatively unchanged untU 
decreases below 750V when the generated plasma becomes unsustainable and 
eventually extinguishes. This is iUustrated m Figure 5 v^^ere fee normalized electron 
20 density shown by diamonds 9 and flie normalized plasma power density shown by 
circles 10 are plotted against the normahzed magnitude of the applied voltage, V JV». 
The normalizations of electron density and plasma power density in this figure, and in 
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Figures 6 and 7 are carried out with respect to the electron density and plasma power 
density obtained by sinusoidal excitation, as iUustrated in Figure 2). Lines 11 and 12 
are drawn through the electron density data 9 and the plasma power desnsity data 10 
respectively, to guide the eye. It is evident that when Vp/Vs > 0.5 plasma pulsing does 
5 not significantly affect electron density while the plasma power is reduced by up to 
40%. Furflier numerical simulations suggest fliat ion densities are also largely 
unaffected as long as Vp/Vs > 0.5. Given liiat electron and ion densities remain 
approximately the same and that the dischai^e current and the gas voltage undergo 
relatively small changes, the basic plasma characteristics should remain relatively 
10 unchanged. In other words voltage pulsing in Figure 4 is unlikely to affect adversely 
plasma-activated applications that rely on the production of electrons oE^and ions, yet is 
capable of considerable power saving. 

It is found that the amount of power saving is dependent on that of nitrogen 
impuiies. In the case of very Uttle nitrogen impurities (<0.01%), numerical studies 
15 suggest a power saving of more than 50% as indicated in Figure 6. Again such 
significant power saving is achieved without affecting the production of electrons and 
ions. Also found is that the power saving is not at the expense of the production of 
. excited neutral species either. With 0.5% nitrogen. Figure 7 shows the normalized 
density of He(2^S) 13 as a function of the Vf/Vg ratio. It is evident that the trend of the 
20 metastable density 13 effectively tracks that of ftie electron density 14. So if the plasma 
power is reduced with the electron density 14 kept approximately the same as that with 
sinusoidal excitation (without pulsing), the metastable density 13 will also remain at 
approximately flie same level as its value wifli sinusoidal excitation. 

The density tracking between electrons and metastables suggests a key role 
25 played by elections in mfluendng production of excited neutral species in diffuse 



atmospheric gas discharges. It is known that energy required to excite neutral species 
and create chemically, reactive species is usually provided through the kinetic energy of 
electrons. There are also reactive species whose densities are directly proportional to 
electron density. For example, densities of most oxygen species (e.g. atomic oxygen, 
5. singlet-sigma metastable oxygen, singlet-delta metastable oxygen) in a diffuse 
nonthermal atmospheric discharge in an oxygen-helium mixture are foxmd to 
proportional to where /» = 0.8 - 1.9 and Pjj^is the input RF power to the generated 
plasma (see J. Y. Jeong, J. Park, I. Renins, S. E. Babayac^ V. J. Tu, G. S. Selwyn, G. 
Ding, and R, F. Hicks, "Reaction chemistry in the afterglow of an oxygen-helium 
10 atmospheric plasma'', J. Phys. Chem., vol. 104, pp. 8027 - 8032; 2000). According to 
eq.(3), these oxygen densities are also proportional to electron density. Therefore the 
very similar electron densities in the sinusoidal case of Figure 2 and the pulsed case of 
Figure 4 may be a good basis to consider that the densities of most reactive species in 
the two cases would be similar and so would their chemical reactivity. Thus it is 
15 possible to produce at least some reactive species electrically efficiently with voltage 
pulsing, and as such this will also be useful for applications that reply on reactive 
species. It is important to note however that.in general there is no direct correlation 
between electron density and densities of all reactive species. There are reactive species 
whose densities may be inversely proportional to electron density, for example ozone in 
20 the oxygen-helium plasma mentioned Jeong et al (supra). So for individual dififuse 
nonthermal atmospheric plasmas, the suggested maintenance of production of reactive 
species deduced from tiaat of electron production needs to be assessed carefully. Our 
observation is that it is possible for diffuse nonthermal atmospheric plasmas to Tnaintam 
densities of some reactive species with reduced electric power. 
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Example Z-^Tail-THmmed Waveform 

An alternative embodiment to achieve the benefits of the iavention may be 
referred to as a tail-trimmed pulse waveform. With reference to Figure 8, such a 
waveform may be described by the following equations: 



F„sin(|-ac,)siii(2«:/'^+^) 



^^+^1 -x^<rem(2rr f t, 27c)<fc, 
{^-x^<remi2}r f t. 2a:) -x, j, 

l^^+x, -Xj J < rem{2 n f t, 2fi) < 2k 



elsewhere 



where : 



\2 TV TV J 



e^ceil{2ft)7t 



and where c«i7(2 / 1) rounds (2 / ^) to the nearest integer towards infinity. 

With reference to Figure 9a, the "tail-off* phase of the applied voltage in 
(between 2S^s and SO^s) may be trimmed to further enhance energy efficiency for 
plasma genemtion. The voltage-current characteristics of the induced atmospheric 
10 plasma from this waveform are shown in Figure 9b, where 15 is the applied voltage, 16 
is the gas voltage, and 17 is the discharge current It is evident that moderate voltage 
t rimming does not significantly affect the generation and characteristics of induced 
nonthermal atmospheric plasmas, though the discharge current 17 has different peak 
values in different half cycles. Numerical studies suggest that further power saving is 
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possible but very modest, typically a few percent and at best 10% before electron 
density starts to decrease. 

We have so far established fliat generation of diffiise nonthennal atmospheric 
plasmas can be made more efficient by altering the waveform of the applied voltage 
without affecting the basic plasma chamcteristics. 

Example 3-Tafl-Shaped Waveform 

In this example, we consider how pulse width may affect plasma generation. la 
this example we consider one particular type of pulsed voltage signal constructed from 
a sinusoidal signal for flie voltage rise phase and a Gaussian decay for the voltage tail 
phase. Matixematically in each cycle tiiey may be expressed as follows 



F(0 = 



V^smot if 2Nn<t<t^^2N7t 

VpQ^^{t-t^^2N7tf lr''\ if t^+2N7r<t<2N7V'¥7C 

V^sraot if 2N7v + 7v<.t<t^+2N7V + 7V 

-F^exp[-(/-^,-2iy^-^)Vr^] if U^2N7v^7t <.t<2Nn^-27t 



where Po is the peak voltage of the sinusoidal signal, is tiie peak voltage of the 
Gaussian decay signal^ t is the pulse width of the Gaussian signal, g> = 2TdT is the 
angular frequency of the applied voltage with T being its period, and /b the instant at 
which the sinusoidal and the Gaussian signals joint. To and U are defined as follows: 



T 

= — arcsin 
2;r 



?o=ro+fcr/2 (5) 

To compaife to the sinusoidal case of Figure 2, we set Vq = 1.5kV, Vp = 1.4kV, and 
T=r/8. The waveform of the applied voltage 18 is shown in Figure 10a, and the voltage 
and current signals of the generated atmospheric plasma is shown in Figure 10b, where 
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19 is liie appUed voltage. 20 is the gas voltage, and 21 is tbe discharge cun«nt It is 
evident that the pattern of one discharge every half cycle remains and this is repetitive 
for many tens of cycles. On the o&er hand, the gas voltage is markedly different from 
that in the sinusoidal case of Figure 2. Particularly it has a large hump preceding the 
6 main peak that causes the discharge event in each half cycle. Also after Ihe discharge 
event where the discharge current peaks, the gas discharge does not reduce to zero as 
rapidly as in the sinusoidal case. 

The discharge current 21 is between 30 - 45mA, markedly lower than 55mA 
calculated for the sinusoidal case. Mterestingly the pulse width of the discharge current 
10 21 is much larger than that in the sinusoidal case (Figure 2), particularly for the positive 
half cycles. For the case in Figure 10. electric power density consumed in the plasma is 
about 0.32WW, about 7% above that in the sinusoidal case. Further numerical 
examples studied suggest that under other pulsing conditions diffiise atmospheric 
heUum plasmas excited with the pulsed voltage of Figure 10a consume approximately 
15 the same amount of electric power as those with sinusoidal excitation. Therefore 
shortening the pulse width may not be as effective in achieving energy saving on an 
absolute basis, but, as discussed below, have advantages in energy-saving for a given 
density of elecUons and metastables achieved. From the voltage-current characteristic, 
tiiis is likely due to the wider pulse widths for both the discharge current and the gas 
20 voltage. On the other hand, electron density achievable can be much higher. For tiie 
case shown in Figure 10b, the peak electron density is found to be 1.74xlO'W^ 83% 
above 0.95x10? W for the sinusoidal case. The spatial variation of electron density 
over one half cycle shows a similar distribution to that that produced by the sinusoidal 
waveform presented above yet with clear difference near the voltage polarity reversal 
25 point. There is a clear evidence of pulsed production of electrons. Further nmnerical 
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calculations suggest that by changing the pulse width of the applied voltage in Figure 
10a, electron and metastable densities can be made much greater than that under 
sinusoidal excitation. Though inappropriate wave-shaping of the applied voltage can 
also reduce electron and metastable densities, flie above discussions suggest that 
5 through control of pulse widfli it is possible to enhance production of electrons and 
metastables for a given plasma pow^ consumption. To see Ms more clearly, we 
calculate the ratio of electron density to plasma power density (number of electrons that 
can be produced at a given input power). This is 5.4xlO^°/W for this example, and 
3.2xlO^^/W for the sinusoidal waveform presented above, an improvement of 68%. 
10 From this comparison, it is still electrically efiScient to pulse the applied voltage for the 
generation of diffuse nonthermal atmospheric plasmas. 

Example 4 - Controlled Waveform 

It is clear that an optimal saving of energy in the production of gas plasmas may be 
15 obtained by the use of a control system, such as a feedback control system, wherein the 
waveform of the applied voltage is driven by the control system, which itself uses the 
discharge current as an input variable. In this way, the voltage signal may be attenuated 
following the plasma discharge, so as to make most efficient use of the input energy. 
The use of such a control system can then also take account of any temporal changes in 
20 the system configuration (such as changes in electrode or dielectric coating 
characteristics, or changes in the gas composition), leading to waveforms (as defined 
herein) that may vary in their repetition rate, and which may vary in peak shape and 
amplitude j&om cycle to cycle. 
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Parametric Ranges of Key System Variables 

Electrode size: Typically between 50cm^ and lOOcm^ of either circular or rectangular 
shapes. Scaling up is possible, with the electrode area divided into different sections 
connected in parallel to the supply power source. To control plasma stability with large 
5 electrodes, it may be useful to use individually valued series resistors or impedance 
networks in each of these parallel circuit branches. On the other hand, plasma stability 
is easier to achieve with smaller electrodes than 50cm^. 

Electrode separation: Typically around 1cm. Smaller electrode separation makes it 
easier to control plasma stability. Although larger electrode separation tends to un- 
10 stabilise the generated nonthermal plasma, it is quite possible to increase the electrode 
separation up to tens of centimetres, especially when the electrode area is large and the 
electric power pulsed. 

Capacitance of dielectric coatings: For the case where at least one of the electrodes has 
a dielectric coating, their capacitance depends on the size of the electrodes and hence 

15 tiiat of the dielectric. Typically for a surface area of lOcm^ the capacitance of one 
dielectric coating is between lOpF and .lOOpF. If the surface area of the dielectric 
coatings is increased, the capacitance increases proportionally. Also it is equally 
acceptable to coat either one electrode or both electrodes, the latter of which will have 
the coating capacitance halved. 

20 Peak voltage: This depends on the repetition frequency, the electrode separation, and 
gas composition. At a repetition fiequency of lOkHz and an electrode separation of 
0.5cm, the peak voltage needs to be greater than l.OkV but usually less than 2kV in 
helium or helium mixed with a small ftaction of nitrogen or/and oxygen (typically less 
than 1%). This needs to be increased when the electrode separation is increased. At a 

25 greater repetition frequency on the other hand, the peak voltage required becomes 
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smaUer. For example between 10 and 20MHz, the peak voltage can be reduced to as 
low as 500V in helium with 0.5cm electrode separation. 

Gas composiHon: Our numerical work was perfoimed for heUum-nitrogen mixture. For 
applications, heUum-oxygen is more useful and argon-oxygen is even better as argon is 
5 cheaper than helium. The most practically convenient gas is air (though very difficult to 
simulate). Alternatively noble gases mixed with a small fraction of a reactive gas (such 
as, but not limited to, oxygen can be used. 

to in the Gaussian decay Junction: to is given in equation 4 and 5 in Example 3 of the 
prefened embodiments. Typically it is between 5% and 15% of the repetition period of 
10 the applied electric voltage. 

Smtunary 

Difiuse nonliiennal gas discharges generated at atmospheric pressure have found 
increasing appKcations in many key materials processing areas such as etching, 
15 deposition, and structural modification of polymeric surfaces. To facilitate tailored 
and improved appKcations of these novel gas plasmas, we consider their pulsed 
generation based on one-dimensional numerical simulation of heUum-nitrogen 
discharges. We consider four waveforms of the plasma-generating voltage, namely 
(1) sinusoidal; (2) peak-levefled sinusoidal; (3) peak-levelled and tail-trimmed 

20 sinusoidal; and (4) pulsed with a Gaussian-sh^ed tail, aU at the same repetition 
frequency of 1 OkHz. For each case, voltage and current characteristics are calculated 
and flien used to assess whether the generated plasma is diffuse and nonthermal. 
Densities of electxMis, ions, and metastables are calculated, together with tiie 
dissipated electric power in the plasma bulk. It is found that plasma pulsmg can 

25 significantly reduce the electiic power needed to sustam difEuse nontiiermal 
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atmospheric plasmas. SpecificaUy by choosing appropriate pulse shape, the plasma- 
sustaining power can be reduced by more flian 50% without reducing densities of 
electrons, ions, and metastables. On the other hand, electron density can be 
enhanced by 68% with the same input electric power if the pulse width is suitably 
5 narrowed. 

Scope of the Inventioii 

TTie invention is defined in the Claims, which now follow, and in which the term 
•Vavefoim" is understood to inchide periodic signals that have a zero value over 
10 more than an instantaneous part of the period, and which are more commonly 
referred to as "pulse trains". 

In the aaims, the term "truncated" is taken to include limited in amplitude to a 
maximum and/or minimum voltage, such limitation typically leading to a waveform 
with a substantially flat profile at its extreme value or values, and/or to include 
15 .narrowed in pulse width. 

la. the Claims, Equation El is taken to be: 
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where: 



and: 



F„8in(|-;^) 



' 2 



elsewhere 



rem{2 n f t, 27v) gives tiie remainder after &e division (2 / )/2;r . 
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and furOier, Equation E2 is taken to be: 



Vif): 



F„sm(|-x,) 

y« sin(j- :r,)sm(2 ?r /■ / + 6*) 
K sin(|~ j:,)sm(2 jrf't + 0) 



K„sm(2«r/0 



0 ^2 ^ 2jc, 

^ TF J 



e = ceil(2f 



^^'^Xy-x^<remi27u f t, 2jc)<2jt 



elsewhere 



nd where cei7(2 / 1) rounds (2 / r) to the nearest integer towards infmity. 
5 and Equation E3 is taken to be: 



sind?^ if 2N7V ^t<tQ+ 2N^ 

V^ejqp[-(t-t,-2N7i:y /r^] if t^^2Nn <.t <21^n ^-n 

V^^m.a)t if 2N7V'{'7C<.t<t^+2N7V + 7V 

^F^exp[-(^-/o -2//^-/r)VT*] if t^+2N7r + 7t <2N7v^2^ 



In all of which equations El, E2 and E3, the parameters are to be interpreted 
described herein. 
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Claims- 

1 . A method of generating a gas plasma characterised by Ihe feature that the appKed 
voltage exhibits a waveform which is truncated. 
5 2. A method of generating a gas plasma characterised by the feature that the appUed 
voltage exhibits a waveform which decays asymmetrically from its peak value. 
3. A melhod of generating a gas plasma characterised by the feature that the appUed 
voltage exhibits a waveform which is truncated and which decays asymmetricaUy from 
its peak value. 

10 4. The method of any of the preceding claims wherein the applied voltage, V, as a 

function of tune, t, said time t being measured from any arbitrary instant^ takes the form 
of a waveform, V(t), of cycle time T. wherein in at least one of the half cycles, i.e. 
between (tr=i T) and (1r=i T + T/2) or between (t?=i T f T/2) and (tKi+l)T), where i takes 
integer values, the waveform is characterised by the magnitude of the integral of the 

15 voltage with respect to time being greater in the first half of said half cycle than in the 
second half of said half cycle. 

5. Ihe method of any of Claims 1 to 3 wherein the appHed voltage. V. as a function of 
time, t, said time t being measured from any arbitrary instant takes the form of a 
waveform, V(t), of cycle time T, wherem in at least one of the half cycles, i.e. between 

20 (t=i T) and (1H T + T/2) or between (tH T + T/2) and (lKi+l)T), where i takes integer 
values, the waveform is characterised by a period of substantially constant voltage. 

6. The method of Claim 5 wherein the appUed voltage is defined by equation El herein. 

7. The method of Claim 4 wherein the appKed voltage is defined by equation E2 herein. 

8. The method of Claim 4 wherein the appUed voltage is defined by equation E3 herein. 
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9. The method of any of Claims 1 to 3 wherein the applied voltage is generated by the 
action of a control system, said control system using a measurement of the plasma 
discharge current as an input signal. 

10. A method of generating a gas plasma substantially as described herein with 
reference to and as illustrated in the accompanying drawings. 
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Abstract 



Generation of DiflEiise Non-Thennal Atmospheric Plasnias 

A method of generating non-tbiermal atmospheric gas plasmas using an 
5 applied voltage with specially-shaped waveforms. The waveforms are shaped 
to avoid energy wastage by attenuating the voltage in the period after current 
discharge has occurred. A control system may also be used to shape tiie 
applied voltage waveform in response to the gas discharge current 



10 Figure 10 should &coowpBiny the abstract 
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Figure 4 
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Figure 6 
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Figure 10. 
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